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FOREWORD 

The  work  described  in  this  report  is  part  of 
the  continuing  study  of  hydrodynamics  at  this 
Station.  The  report  analyzes  the  effect  of  body 
shape,  tail  and  rudder  size,  and  other  physical 
parameters  on  the  stability  and  response  of  tor¬ 
pedoes.  The  results  of  this  analysis  will  be 
useful  in  preliminary  design  and  control  studies. 

The  work  was  performed  between  June  1957 
and  June  1958  under  Bureau  of  Ordnance  Task 
Assignment  NO-404-664/4100  l/O  1060.  The  re¬ 
port  was  reviewed  for  technical  accuracy  by 
H.  T.  Yerby  and  L.  A.  Lopes  of  this  Station. 


D.  J.  WILCOX,  Head 
Underwater  Ordnance 
Department 


Released  under 
the  authority  of: 

WM.  B.  McLEAN 
Technical  Director 
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ABSTRACT 

The  stability  and  response  of  torpedoes  are 
affected  by  variations  in  body  shape,  weight, 
volume,  center  of  gravity,  tail  size,  rudder 
size,  and  velocity.  Two  hydrodynamic  design 
criteria  can  be  developed  which  provide  ade¬ 
quate  stability  and  good  response.  All  torpedoes 
of  conventional  shape,  designed  to  the  criterion 
L*,  a  dimensionless  tail-lift  coefficient,  respond 
to  stepwise  rudder  deflections  in  an  almost  iden¬ 
tical  manner,  provided  that  the  steps  result  in 
equal  changes  in  the  steady- state  turn  rate.  To 
obtain  such  almost  identical  response  in  the 
special  case  of  linear  proportional  control,  a 
second  criterion  K*,  characteristic  for  control 
gain,  is  required.  These  two  criteria  make  it 
possible  to  match  the  response  of  widely  different 
torpedo  bodies  both  to  constant- rudder  and  to 
linear  proportional  control. 
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NOMENCLATURE 
Maximum  cross-sectional  area,  ft 

Constants  defined  in  Appendix  A 

■¥ 

Prismatic  coefficient,  c  =  — 

p  At 

Drag 

(1/2)PAV2 

( l/2)pAf2  J^,  moment  of  inertia  about  z  axis,  slug  ft2 
A  constant  defined  in  the  control  equation 
A  particular  value  of  K 

A  particular  value  of  defined  on  page  3 

Tail-lift  derivative 
(l/2)pAV2 
Total  length,  ft 
Cylindrical- section  length,  ft 
Nose  length,  ft 

*n  +  lt  >  ft 
Tail  length 

mj  and  -go£,  ±  w V  1  -  £  ^  i  are  the  roots  of  the  cubic 
equation  associated  with  the  yaw  equation  with  linear 
proportional  control 
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6r 

Nr 
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YP 


A  particular  value  of  m] 

Apparent  longitudinal  mass  (Ref.  2,  pp.  22  -  23) 

IpAfm^ 

Apparent  transverse  mass 
■|pAf  m'rj. 

Moment  derivative  with  respect  to  p,  ft  lb 
( l/2)pAiV2N  pg  =  body- moment  derivative,  ft  lb 

Rudder-force  derivative  xg  ( 

(i/2)PAiv2  e  6r 

Damping-moment  derivative,  ft  lb 

( l/2)pAf  2VN'rg  =  body- damping-moment  derivative, 
ft  lb 

z- component  angular  velocity 
Time 

Speed,  ft/sec 
Volume  of  torp  >do,  ftJ 

Right-hand  coordinate  system,  fixed  in  torpedo  body, 
originating  on  torpedo  axis  at  the  point  nearest  the 
center  of  gravity;  x-axis  along  torpedo  axis  positive 
forward,  z  positive  downward 

Fixed,  right-hand  coordinate  system  with  zQ  positive 
downward 

Distance  between  origin  of  moving  coordinate  system 
and  center  of  rudder  lift,  ft 

Distance  between  origin  of  moving  coordinate  system 
and  center  of  tail  lift,  ft 

Force  derivative  with  respect  to  p,  lb 
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yPb 


f 

(l/2)pAV*Y  Pg  =  body-force  derivative  with  respect 
to  p,  lb 

Rudder-force  derivative 
(l/2)pAV2 

Damping-force  derivative,  lb  sec 

( l/2)pAf  VYrg  =  body-damping-force  derivative, 
lb  sec 


P  Angle  of  attack  in  yaw  =  angle  between  the  velocity 
vector  and  its  projection  on  the  x-z  plane  (measured 
in  direction  from  the  projection  toward  the  y  axis) 

6r  Rudder-deflection  angle,  positive  in  direction  from 
-x  axis  to  -y  axis,  rad 

£  Damping  factor 

■H  ^n/  -*p 

0  Pitch  angle:  angle  between  projection  of  x  axis  on 

x0  -  yQ  plane  and  x  axis  (measured  from  xD  -  yD  plane 
in  direction  of  -zQ  axis),  rad 

p  Density  of  fluid,  slugs/ ft2;  p  ~  2  slugs/ft^  for  salt 
water 

<|>  Roll  angle:  angle  between  plane  containing  x  and  zQ 
axes  and  -z  axis  (measured  from  the  plane  to  the 
-z  axis  in  direction  from  -z  to  y),  rad 

iJj  Yaw  angle.-:  angle  between  xQ  axis  and  the  projection 
of  x  axis  on  xQ  -  yQ  plane  (measured  from  xQ  in 
direction  of  yQ),  rad 

dijj 

^  — 
dt 
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INTRODUCTION 

The  ability  of  a  torpedo  to  pursue  a  target  successfully  re¬ 
quires  that  an  error  in  heading  be  corrected  quickly,  without 
objectionable  overshoot  or  oscillation;  these  corrections  should 
be  made  by  controls  which  are  not  excessively  large  or  powerful. 
This  report  presents  hydrodynamic  design  criteria  for  a  wide 
variety  of  body  shapes  which  give  adequate  stability  and  have 
good  response  characteristics.  Moreover,  the  response  of  all 
torpedoes  designed  to  the  criteria  will  be  almost  identical,  so 
that  the  question  of  matching  the  response  of  widely  different 
bodies  is  at  least  partially  answered. 

A  group  of  ten  torpedo  bodies,  with  the  most  important  body 
parameters  varied  systematically  to  include  most  practical  con¬ 
figurations,  was  studied  to  show  the  effect  of  body  shape  on 
stability  and  response.  The  investigation  was  limited  to  the  yaw 
plane,  but  the  response  in  the  pitch  plane  is  similar,  except  for 
the  effects  of  weight  and  buoyancy. 


CONTROL  SYSTEM 

In  order  to  develop  the  desired  design  criteria,  a  type  of 
control  system  has  to  be  selected.  One  of  the  simplest,  insofar 
as  instrumentation  and  circuitry  are  concerned,  is  linear  pro¬ 
portional  control.  If  vp  is  the  yaw  error,  the  rudder  deflection 
6r  is  given  by  6r  =  Ci(j,  where  C  is  a  constant  called  the  control 
gain  (see  Nomenclature).  Analysis  of  a  linear  proportional  con¬ 
trol  system  may  lead  to  excessive  component  sizes  and  suggest 
a  non-linear  syste  n  or  the  incorporation  of  derivative  control. 
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However,  because  derivative  control  or  feedback  from  control 
surfaces  introduces  added  instrumentation,  networks,  and  ampli¬ 
fiers,  the  use  of  linear  proportional  control  was  assumed.  This 
assumption  may  not  be  too  restrictive,  since  a  torpedo  having 
good  control  characteristics  with  this  system  may  be  expected  to 
respond  adequately  to  other  types  of  control.  For  the  torpedo 
bodies  considered  in  this  report,  linear  proportional  control 
results  in  reasonable  hydrodynamic  design  requirements. 


DESIGN  CRITERIA  FOR 
ADEQUATE  STABILITY  AND  GOOD  RESPONSE 


According  to  Appendix  A,  if  the  control  equation  is 

K 

(1)  6r  = - >|; 


then  the  yaw- error  equation  is 

(2)  i  +  (a  +  bLpT)i  +  (c  +  dL’pT  +eK)i  +  (Kf)iji  =  0 

where  a,  b,  c,  d,  e,  and  f  are  constants  for  a  particular  torpedo, 
as  defined  in  Appendix  A.  Equation  2  has  the  general  solution 

(3)  4*  =  A  exp(mjt)  +  [B  cos  (uV*  -  ^  t) 

I  y 

+  C  sin(co\l  -  £“t)]  exp(-u)t,t) 

where  mj  and  ±  coVl  -  £^i  are  the  roots  of  the  auxiliary 
equation  associated  with  Eq.  2.  Since  ip  is  to  approach  zero,  rrq 
must  be  negative,  for  stability.  The  quantities  u  and  are 
called,  respectively,  the  natural  frequency  and  the  damping 
factor.  If  w  and  £  are  specified,  then  (from  Appendix  B) 
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2 e£u2  +  (f  -  4f £2  -  aeju1”  +  2af£,o  -  cf 

<*».  0  =  - 7 - 

bea)^  -  2bft,w  +  df 

bw^  -  2d£,w2  +  (da  -  bc)w2 


(5) 

k(«,  t,)  = - - - 

bet/  -  2bf£w  +  df 

and 

(6) 

-fK(co,  i) 

m/w.  i)  =  7 

01 

where 

tail  lift 

PT  {l/2)pAV2 


where  p  =  density  of  water.  Appendix  D  gives  the  derivation  of 
Lp.j.  as  a  function  of  fin  size.  The  values  w  =  2tt  and  £  =  0.707 
are  considered  quite  satisfactory  for  frequency  and  damp¬ 
ing.  *  The  corresponding  dimensionless  tail-lift  coefficient, 
L*=  Lprp(2ir,  0.707),  and  the  constant  K  from  the  control  equa¬ 
tion,  K*  =  K(2ir,  0.707),  will  then  be  taken  as  the  design  criteria. 
If  also  m*  =  m^(2ir,  0.707)  then,  from  Eq.  4,  5,  and  6 

f(8.886a  -  c  -  39.48)  +  e(350.7  -  39.48a) 

(7)  L*  = - 

f(d  -  8.886b)  +  39.48be 
b(  1558  -  39.48c)  +  d(39.48a  -  350.7) 


(8) 

K*  = - 

f(d  -  8.886b)  +  39.48be 

and 

(9) 

-fK* 

mf  - 

*  o  r\  A 

39.48 


*  These  values  were  suggested  by  L.  A.  Lopes  and  D.  E.  Elliott 
on  the  basis  of  past  experience. 
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For  these  values  Eq,  3  becomes 


/  -fK*  \ 

( 10)  i}i  =  A  exp  I - t 

\ 39.48  / 

+  (B  cos4.443t  +  C  sin4.443t)  exp(-4.443t) 


If,  as  is  often  the  case,  the  rudder  deflection  is  limited  to 
the  range  {-5Vq,  6 rQ )  then,  from  Eq.  1,  the  range  in  4*  (referred 
to  as  the  "band  of  proportional  control")  for  which  Eq.  9  applies 
is 


(ID 


Thus  the  criteria  L*  and  K*  lead  to  good  control  characteristics 
for  small  yaw  errors.  If,  however,  the  magnitude  of  the  yaw 
error  i|j  is  greater  than 


r 


o 


K* 

the  rudder  will  be  "hard  over"  and  the  resulting  yaw  response 
will  be  of  importance.  The  yaw  equation  is  then 


(12)  iji  +  (a  +  bL*)(p  +(c  +dL*)ijj  = 


The  steady-state  solution  of  this  equation  is 
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(13) 


* 


ss 


6 

r 


rQ 


=  + 


c  +  dL* 


Here  the  main  questions  are  whether  or  not,  for  a  given  torpedo 
body  and  a  given  steady- state  turn  rate,  the  design  criterion  L* 
will  lead  to  a  reasonably  quick  adjustment  to  the  steady- state 
turn  rate,  and  how  the  response  times  will  differ  for  two  different 
bodies.  It  will  be  seen  later  that  the  answers  to  both  questions 
are  favorable;  these  criteria  give  good  results  for  large  errors 
(hard-over  rudder)  as  well  as  small  ones  (proportional  control). 


PARAMETERS  OF  SAMPLE  TORPEDOES 

A  group  of  ten  torpedo  bodies  is  defined  in  Table  1  and  shown 
in  Fig.  1.  The  nose  of  each  torpedo  is  assumed  to  be  elliptical 
and  the  tail  parabolic,  according  to  equations  given  in  an  earlier 
report.^  If  a  nose  has  a  flat  portion,  its  effect  is  assumed  to  be 
too  small  to  be  significant.  The  three  parameters  r)  =  *n/(£n  +  it), 
ip/d  =  (in  +  it)/d,  and  ic/i,  defined  in  Table  1,  are  then  suf¬ 
ficient  to  describe  the  torpedo  shape.  These  parameters  are 
varied  independently  through  ranges  which  include  most  prac¬ 
tical  torpedo  configurations,  so  that  conclusions  concerning 
these  bodies  will  be  generally  applicable.  Additional  assump¬ 
tions  are: 

1.  Speed  V  =  40  knots. 

2.  Volume  if  =  5  ft3. 

3.  Specific  gravity  1.2,  weight  uniformly  distributed. 

2 

U.  S.  Naval  Ordnance  Test  Station.  A  Fast  Method  for 
Finding  the  Physical  Characteristics  of  a  'Torpedo,  by  D.  Argue. 
Pasadena,  Calif.  ,  22  August  1956.  (G.  C.  Memorandum  685.) 
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TABLE  1.  Physical  Parameters  of 
Sample  Torpedoes 


Body 

•^n 

ip  in  4  it 

•^c 

i 

d 

i 

q  - 

in  +  h 

d  d 

i 

1 

0.1 

3 

0 

3 

4.716 

1.1 

0.1 

3 

0.6 

7.5 

7.590 

2 

0.1 

5 

0 

5 

6.630 

2.1 

0.1 

5 

0.6 

12.5 

10.675 

3 

0.4 

3 

0 

3 

4.605 

3.1 

0.4 

3 

0.5 

6 

6.612 

4 

0.4 

4 

0 

4 

5.580 

4.1 

0.4 

4 

0.7 

13.3 

10.893 

5 

0.4 

5 

0 

5 

6.475 

5.1 

0.4 

5 

0.5 

10 

9.290 

d  =  diameter 

f  =  total  length 

lc  =  length  of  cylindrical 
section 

fn  =  length  of  nose  section 


it  =  length  of  tail 
section 

^p  =  ^n  *  ^t 

V  =  40  knots 
-¥  =  5  ft3 


4.  The  dimensionless  drag  coefficient,  based  on  maximum 
cross-sectional  area,  is  0.1  (this  is  not  a  critical  quantity). 

5.  xg  =  1.05  Xpj  Xp  is  the  distance,  in  feet,  from  the  CG  to 
the  center  of  tail  lift  (xp  <  0);  xg  is  the  distance  from  the 
CG  to  the  rudder  center  of  lift. 

6.  The  tail  lift  acts  at  a  distance  i  -  0.293ij.  -  0.167  feet 
from  the  nose.  This  is  approximately  the  point 
where  the  body  diameter  is  one  half  the  maximum 
body  diameter. 

7.  When  the  rudders  are  deflected,  there  is  a  certain  time 
lag  before  the  full  lift  develops.  According  to  unsteady- 
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airfoil  theory,  this  time  lag  approximates  the  time  required  to 
travel  three  or  four  chord  lengths  or,  for  all  cases  considered 
here,  about  0.01  -  0.02  second.  Therefore  this  effect  was 
neglected. 

The  physical  characteristics  and  dimensionless  hydrodynamic 
body  coefficients,  necessary  for  the  determination  of  the  quan¬ 
tities  a,  b,  c,  d,  e,  and  f,  are  found  for  each  body  from  Ref.  1 
and  2  and  the  report  referred  to  in  footnote  2.  Then  L*,  K*, 
and  m*  are  computed  from  Eq.  7,  8,  and  9.  Table  2  gives 
these  values,  the  presently  used  stability  criterion  G  (Ref.  2) 
and  -Y g r/(4,ss/^r0)  from  Eq.  13.  Appendix  C  presents  the  com¬ 
putation  of  the  criteria  for  one  of  the  bodies  (Body  1.1)  in  order 
to  demonstrate  the  method  by  a  practical  example.  It  turns  out 
that  the  smallest  value  of  -m’J  is  11.1,  so  that  in  each  case  the 
effect  of  the  root  m*  can  be  neglected  after  a  short  period  of 
time  (see  Appendix  B).  Other  conclusions  are:  For  a  given 
nose  and  tail  shape,  both  L*  and  -K*  increase  by  adding  a 
cylindrical  section.  The  addition  of  a  cylindrical  section  in¬ 
creases  the  Ygr  required  for  a  given  steady- state  turn  rate. 

The  effect  of  the  relative  lengths  of  the  nose  and  tail  sections  may 
be  assessed  by  comparing  Bodies  1  and  3  or  2  and  5.  As  the 
ratio  fnAt  increases  [holding  (fn  +  ft  )/d  constant],  L*  increases 
while  K*  and  Y 5 r/(ips s/6 rQ)  remain  essentially  unchanged.  Con¬ 
versely,  with  a  constant  in/ h>  the  effect  of  increasing  (fn  +  ft)/d 
is  to  increase  L*,  K*,  and  Y^_  / (‘I'ss/ ^  xQ) • 
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CD  CG  Ypfl*07  Y  rjj  »  0  065  m  L  «  0  5464  (2  »p  gr  4  0  26) 

0.3427.  cp.  0.5464.  „  0.606.  N'rB‘-0  03.  m'T  •  0  5464  (2  «p.  gr  4  16).  j',  ■  0  0275  4  0  0598  ip.  gr 


CD  CG  Ya„«0  5 

0.426.  cp  *  0.8186.  ^  .  ,  , 


0  54  Y'rB  *  0  30  m'L  *  0  8186  (2  *p.  gr  4  0.05) 


-0  15.  m'T  •  0.8186  (2  ip.gr.  4  1  82).  j',  ■  0  0815  4  0.0988  ip.  gr. 


CB  CG 

- - - .  0,343.  c_«  0.5464. 


YpB  •  0.73  Y*rB  »  0.063  ■  0.5464  (2  »p,gr  4  0  11) 

NpB  «  0.608,  N'ry  °  -0.03.  m'j  ■  0.5464  (2  «p.  gr.  4  1.78),  j'z  ■  0  0346  4  0  0494  «p.  gr 


5£.^2.0  427  c  .  0  6)86  YfB  "  075  Y'rB  .  0.30  m'L  .  0  8186  )2  .p.gr  4  0  02) 

,  "  ,  ■  P  '  '  N()d.1.I,  N'rD  .  -0  15.  m'T  .  0  8186  (2  ip.gr  4  )  96).  j',  ■  0.0886  4  0  0973  ip  gr 


CB  CG  Y0.  -0.66  Y*rn  "  0  085  m'.  ■  0,5866  (2  >p  gr  4  0  26) 

- - - .  0  434,  c„  •  0.5866,  ,PB  ,  B  ,*>  , 

/  f  P  NpB»077,NrB«-004,mj*  0.5866  (2  «p.  gr.  4  16),  J  z  »  0.0257  4  0  0558  tp  gr 

FIG.  1.  Sample  Torpedo  Bodies. 
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YAW-ERROR  CORRECTION  FOR 
SAMPLE  TORPEDOES 

Figure  2  shows  the  correction  of  a  20-degree  yaw  error  for 
each  of  the  ten  bodies  described  in  the  preceding  section,  under 
the  assumption  that  the  rudder  deflection  is  limited  to  +6  degrees 
and  that  the  steady-state  turn  rate  for  this  deflection  is  +20  deg/sec. 
The  initial  conditions  are 

i)j  =  20  degrees 

4*  =  o 

V2  N's 

ijj  =  6 - deg/sec4 

iZ  J'z 

at  t  =  0.  The  error  is  corrected  in  slightly  more  than  one 
second,  with  very  slight  overshoot.  That  the  design  criteria  L* 
and  K*  lead  to  torpedoes  with  quite  similar  responses  is  readily 
apparent  since  the  curves  es sentially  overlap. 

The  curves  in  Fig.  2  are  plotted  without  regard  for  the  time 
lag  in  the  control,  under  the  assumption  that  the  rudder  is  in¬ 
stantaneously  at  the  angle  prescribed  by  the  control  Eq.  1. 
Practically,  however,  the  rudder  "lags"  behind  the  desired 
angle.  Therefore  the  control  is  more  accurately  described  by 
Tp6r  +  Sr  =  K*/(.Xg/f) where  Tp  is  the  time  lag.  Figures  3 
through  12  ,  show,  for  Lp^  =  L*  and  K  =  K*,  the  effects 
of  a  0.05-second  and  0.25-second  time  lag  on  each  body.  The 
0.05-second  lag  produces  no  great  effect  on  the  response,  re¬ 
sulting  only  in  slightly  more  overshoot.  A  time  lag  of  0.25  sec¬ 
ond,  however,  results  in  relatively  wide  and  prolonged  oscilla¬ 
tions.  Therefore  the  time  lag  should  be  held  to  a  minimum,  on 
the  order  of  0.05  second  or  less. 

Thus  far  it  has  been  demonstrated  that,  independent  of  body 
shape,  the  criteria  L*  and  K*  result  in  torpedoes  with  very 
similar  responses,  under  the  assumptions  stated  in  the  preceding 

1  1 
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FIG.  2.  Correction  of  Yaw  Error  for  Each  of  the  Sample  Bodies 
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IG.  7.  Response  of  Body  3,  Including  Effect  of  Time  Lag  (Tp)  and  Changes 
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section.  The  question  arises  as  to  how  changes  in  the  assump¬ 
tions  would  affect  the  criteria  and  the  response.  Table  3  shows 
the  results  of  independent  changes  of  Body  1.1  in  weight,  volume, 
speed,  and  location  of  the  center  of  tail  lift  (i,  e.  ,  Xp /t).  The 
changes  are  specified  in  the  table  and  values  L*,  K*,  and 
~Y  Sj-Z^ssArQ)  are  given.  Figure  13  shows  the  correction  of  a 
20-degree  yaw  error  for  each  of  these  cases;  the  obvious  con¬ 
clusion  is  that  the  response  is  not  affected  as  long  as  the  dimen¬ 
sionless  tail-lift  coefficient  and  the  constant  K  from  the  control 
equation  are  given  by  the  criteria  Lp^,  =  L*  and  K  =  K*  (assuming 
the  same  steady- state  turn  rate). 

It  might  also  be  important  to  discover  the  effect  of  replacing 

3 

the  assumption  of  a  5  ft  volume  by  the  assumption  of  a  constant 
cross-sectional  area.  Table  4  shows  that  the  addition  of  a 
cylindrical  section  requires  a  larger  L.*,  -K*,  and  -Y  5 ^/(dj s s/ S r  ). 
The  same  is  true,  although  to  a  lesser  degree,  if  the  point  of 
maximum  diameter  is  moved  rearward,  without  addition  of  a 
cylindrical  section. 

EFFECT  OF  CRITERIA  CHANGES  ON  RESPONSE 

The  effect  of  separate  changes  in  Lp,p  and  K  on  the  response 
of  each  body  is  also  investigated  in  Fig.  3  through  12.  For  the 
recommended  Lp,p  =  L *  and  a  time  lag  of  0.05  second,  K  was 
given  the  two  values  1.5K*  and  O.SK*.  The  response  for  each 
of  these  values  does  not  vary  significantly  from  that  obtained 
with  K  =  K*.  However,  holding  K  =  K*  fixed  and  varying  Lp<j, 
leads  to  quite  significant  changes.  A  20%  increase  in  LpT  (to 
L'pt  =  1.2L*)  results  in  slower  correction  of  a  yaw  error,  while 
a  20%  decrease  (to  Lp,p  =  0.8L*)  leads  to  large-amplitude  os¬ 
cillations.  The  effects  are  more  pronounced  on  bodies  without 
a  cylindrical  section.  Figure  14  shows  that  at  these  changed 
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TABLE  4.  Effect  of  Change  in  Body  Shape  if 
Cross-Sectional  Area  Is  Kept  Constant 


Body  Description 

L* 

BB 

-K* 

‘I'ss/^rQ 

Body  2 

*n  *c 

- =0.1,  —  =  0 , 

fn  +  ~'t  * 

1.31 

0.0352 

0.079 

*n  +  *t 

-  =5,  i  =  6.63  ft 

d 

Body  2  with  cylindrical 
section  added  so  that 

^c 

—  =  0.6 

SL 

1.94 

0.2181 

0.673 

Same  cross-sectional  area 
as  Body  2 

^n  ^c 

-  =  0.4,  -—  =  0, 

*n  +  *t  l 

! 

j 

i 

1.66 

0.0412 

0.089 

*n  +  h 

- =5,  i  =  6.63  ft 

d 

values  of  L  the  response  cannot  be  improved  significantly  by 
varying  K.  This  does  not  mean,  however,  that  L*  is  an  optimum 
value;  the  explanation  can  be  found  in  Expression  11  and  Eq.  13. 

In  Eq.  13  the  steady- state  turn  rate  is  directly  proportional  to 
Ygr  and  inversely  proportional  to  a  linear,  increasing  function 
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of  Lp^.  Thus,  if  Lp^,  is  reduced  to  0.8L*  and  Yg^  is  not 
changed,  the  steady- state  turn  rate  is  increased,  which  accounts 
for  the  faster  initial  response.  Also,  since  the  band  of  propor¬ 
tional  control  has  not  changed  (Expression  11),  the  turn  rate  is 
sufficiently  high  to  cause  the  torpedo  to  overshoot  and  the  rudders 
to  become  hard-over  at  -6  degrees.  In  fact,  the  overshoot  is  so 
great  that  the  torpedo  has  time  to  build  up  a  sufficiently  large 
turn  rate  in  the  opposite  direction  to  again  overshoot  the  propor¬ 
tional  band,  which  results  in  the  extreme  oscillation.  In  the  case 
of  =  1.2L*,  with  no  change  in  Yg^,  the  turn  rate  is  decreased 

so  that  no  overshoot  occurs,  but  the  time  needed  to  correct  the 
error  increases.  If  Yg_  is  reduced  for  Ep^  =  0.8L*  and  in¬ 
creased  for  Lp^,  =  1.2L*,  so  as  to  maintain  the  original  steady- 
state  turn  rate  (Fig.  3  through  12),  the  initial  part  of  the  curves 
will  be  approximately  the  same,  differing  only  because  of  dif¬ 
ferent  lengths  of  time  required  to  reach  the  steady- state  turn 
rate,  In  the  former  case  (Lp^,  =  0.8L*  and  a  reduced  YgJ  the 
band  of  proportional  control  is  narrowed  and  the  danger  of  os¬ 
cillation  by  overshooting  the  band  persists.  K  may  be  decreased 
to  widen  the  band,  but  then  both  and  K  differ  from  L*  and 

K*  ar.d  the  corresponding  frequency  and  damping  (w  and  t,)  may 
be  undesirable.  In  the  latter  case  (L«„  =  1.2JL*  and  an  increased 
Y g^)  the  band  of  proportional  control  is  widened  so  that  there 
should  be  no  danger  of  overshooting  the  band.  In  fact,  K  may  be 
increased  to  reduce  the  proportional  band  width. 

Figures  15  and  16  indicate  how  w  and  £  change  with  L'p^  and 
K.  The  curves  were  plotted  for  Body  1.1  only,  but  corresponding 
changes  on  the  other  bodies  should  resuit  in  similar  curves.  The 
principal  conclusions  are  that  L'pT  is  largely  independent  of  w 
and  depends  primarily  on  C,  while  K  is  independent  of  £  and  de- 

1 

pends  mainly  on  to.  Also,  Lp,p  increases  with  £  and  -K  increases 
with  to. 
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RESPONSE  OF  SAMPLE  TORPEDOES  TO 
CONSTANT  RUDDER  DEFLECTION 

Figure  17  shows  how  each  of  the  ten  bodies  approaches  a 
20  deg/  sec  steady- state  turn  rate  for  Lp,p  =  L*  and  a  rudder  de¬ 
flection  of  -6  degrees.  The  rudders  were  assumed  to  deflect 
from  0  to  -6  degrees  in  0.05  second  at  a  constant  rate.  The 
curves  do  not  differ  greatly,  further  demonstrating  the  constancy 
of  response  with  respect  to  the  criteria,  and  Body  1.1  gives  the 
central  curve.  Figure  18  shows,  for  Body  1.1,  the  effect  of  a 
±20%  change  in  Lp^,  [i.  e.  ,  Lp,p  =  (1  ±  0.2)L*],  with  a  correspond¬ 
ing  change  in  to  give  the  same  steady- state  turn  rate.  A 
change  of  approximately  5%  in  L*  would  make  all  the  curves  of 
Fig,  17  practically  identical. 

Because  of  its  importance  in  determining  cavitation  resistance, 
the  angle  of  attack  at  the  tail  is  plotted  in  Fig.  19  for  each  of  the 
ten  bodies. 

APPROXIMATE  FORMULAS  FOR  DESIGN  CRITERIA 

The  following  approximate  formulas  for  L*,  K*,  and 
■Y  6r/(llJss/6r0),  in  terms  of  the  body  parameters  t]  =  fn/(^n  +  ft)> 
fc/i,  and  fp/d  =  (in  +  it)/ d,  should  give  sufficiently  accurate, 
values  for  design  purposes: 


(H) 


L*~  (0.909  -  0.37ti)  —  +  (r\  +  1.205) 

l 


I  io  r 

(15)  -K*  a  — r|  -  - 
\  3  3 


/  >  \ 

I*  c) 

2 

«c  / 

10.314—  -  0.0131 

l7j 

0.15—  +  0.103 
d 

7  +  ( 

/P 

d 


*1 

d 


f  n  fp 

+  0.019  -  0.016  +  0.1  —  +  0.082  — 
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and 


(16)  7 


-Y6r 


^Bs/^Tc 


/10  1  1 
*1 


0.052—  +  0.067 
d  / 


—  -  0.0235—  +  0.0512  —+  0.0019 -  0.0051 


+  0.007— +  0.0185-- +  0.0448 
d  d  / 


These  approximations  save  considerable  computing  time  with  quite  satis¬ 
factory  accuracy.  Table  5  compares  the  approximate  values  for  the  ten  bodies 
under  consideration  with  the  exact  values.  Two  other  bodies,  with  parame¬ 
ters  t],  fc/f,  and  fp/d  lying  within  the  range  of  the  ten  sample  bodies,  were 
chosen  at  random.  The  approximate  formulas  hold,  the  largest  error  being 
11%  for  K*. 

These  approximations  were  established  under  the  assumptions  of  the 
section  "Parameters  of  Sample  Torpedoes"  and  should  therefore  be  used 
with  caution  if  a  particular  torpedo  departs  significantly  from  any  of  the  as¬ 
sumptions.  Table  3  shows  the  effect  of  changes  in  certain  of  the  assump¬ 
tions  for  one  of  the  sample  bodies;  if  the  approximate  formulas  are  used 
with  corrections  as  indicated  in  this  table,  good  approximations  should  still 
result. 


CONCLUSIONS 

It  has  been  shown  that  all  torpedoes  of  conventional  form,  designed  to 
the  criteria  L*  and  K*,  will  have  practically  identical  responses  both  for 
the  constant-rudder  case  and  for  linear  proportional  control.  The  intention 
was  to  obtain  a  reasonably  adequate  response  rather  than  an  "optimum" 
response.  In  fact,  the  term  "optimum"  itself  is  questionable;  it  is  not  at  all 
certain  that  a  response  which  corrects  an  error  in  the  shortest  time, 
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regardless  of  the  amount  of  overshoot  and  oscillation,  is  prefer¬ 
able  to  one  which  takes  longer  but  results  in  less  overshoot  and 
oscillation.  This  question  cannot  be  answered  in  general  terms. 

It  can  be  answered  only  for  a  specific  torpedo,  after  physical 
restrictions  and  performance  requirements  are  known  from  com¬ 
puter  and  weapons  analysis  studies.  Past  experience  indicates, 
however,  that  the  response  of  a  torpedo  with  a  dimensionless 
tail-lift  coefficient  L*  is  quite  adequate,  and  certainly  adequate 
for  preliminary  design. 

It  should  be  noted  that  the  quantity  K*  is  uniquely  associated 
with  the  linear  proportional  control  system  and  hence  has  no 
meaning  for  other  control  systems.  However,  consideration  of 
such  a  linear  proportional  control  system  in  this  study  has 
served  the  following  important  purposes:  It  has  lead  to  values  of 
the  dimensionless  tail-lift  coefficient  L*  which  result  in  torpedoes 
with  an  adequate  degree  of  stability  and  approximately  equal  re¬ 
sponse  to  a  constant  rudder  deflection.  Further  it  has  shown 
that  torpedoes  having  Lp,p  =  L*  can  be  controlled  by  a  linear 
proportional  control  system  and  hence  will  probably  behave  well 
under  other  control  systems.  Since  any  motion  of  the  control 
surface  can  be  approximated  by  a  succession  of  constant  rudder 
deflections,  and  since  the  torpedoes  respond  similarly  with  a 
constant  (step)  deflection,  it  should  be  possible  to  make  them 
respond  similarly  with  any  type  of  control. 

Some  of  the  more  important  conclusions  can  be  summarized. 

1.  Bodies  of  Equal  Volume.  An  elongated  streamlined  body, 
or  a  body  with  a  cylindrical  section,  has  a  smaller  steady- state 
turn  rate  per  degree  of  rudder  deflection  than  a  shorter  stream¬ 
lined  body,  if  tail  and  rudder  sizes  are  equal;  to  obtain  equal 
response  from  the  elongated  body  the  size  of  the  tail  must  be 
decreased  and  the  size  of  the  rudder  increased. 
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Moving  the  maximum  diameter  of  a  streamlined  body  rear¬ 
ward  increases  the  steady- state  turn  rate  per  degree  of  rhdder 
deflectioh;  to  obtain  equal  response,  the  size  of  the  tail  must  be 
increased,  but  the  size  of  the  rudder  can  remain  practically  un¬ 
changed. 

2.  Bodies  of  Equal  Maximum  Cross-Sectional  Area.  The 
addition  of  a  cylindrical  section  decreases  the  steady-state  turn 
rate  per  degree  of  rudder  deflection,  if  tail  and  rudder  size  are 
equal;  to  obtain  equal  response,  the  size  of  both  tail  and  rudder 
must  be  increased. 

Moving  the  maximum  diameter  rearward  increases  the 
steady- state  turn  rate  per  degree  of  rudder  deflection;  to  obtain 
equal  response,  the  size  of  both  tail  and  rudder  must  be  in¬ 
creased. 

3.  The  steady- state  turn  rate  per  degree  of  rudder  deflection 
is  a  decreasing  function  of  tail  size  and  an  increasing  function  of 
rudder  size. 

4.  If  the  steady- state  turn  rate  per  degree  of  rudder  deflec¬ 
tion  is  held  constant,  a  torpedo  with  larger  rudder  and  tail  will 
respond  more  quickly  than  a  torpedo  with  smaller  rudder  and 
tail. 

5.  If  the  steady- state  turn  rate  per  degree  of  rudder  deflec¬ 
tion  is  held  constant  by  adjusting  the  size  of  the  rudder,  small 
changes  in  torpedo  weight,  volume,  speed,  and  center  of  gravity 
have  little  effect  on  response,  provided  the  tail  lift  is  changed 
according  to  the  criterion  developed  in  this  study. 

6.  In  general,  a  time  lag  of  0.05  second  in  the  control  system 
appears  to  be  acceptable,  but  a  time  lag  of  0.25  second  apparently 
is  not. 

7.  According  to  Ref.  1  geometrically  similar  bodies  have 
equal  dimensionless  body  coefficients.  If  the  bodies  also 
are  of  uniform  density  with  the  same  specific  gravity  (so  that 
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the  dimensionless  masses  and  moments  of  inertia  are  equal), 
and  if  the  centers  of  tail  and  rudder  lift  are  located  so  that  Xp  / 1 
and  Xg/i  remain  constant,  the  bodies  will  have  identical  re¬ 
sponses  for  both  linear  proportional  control  and  step  rudder  de¬ 
flections  as  long  as  the  ratio  V/i  is  constant.  This  is  due  to  the 
fact  that  the  constants  a,  b,  c,  d,  e,  and  f  (Appendix  A)  are  then 
equal  and  hence  the  equations  of  motion  are  identical. 
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Appendix  A 

EQUATIONS  OF  MOTION  AND  CONTROL 


The  motion  of  a  torpedo  in  yaw  (Ref.  2  or  3)  is  given  by 


Jz4-  -  Nrij>  -  Npp  =  N6r6r 

and 

(VmL  -  Yr)ij<  -  VmTp  -  Ypp  =  Y6r6r 


These  are  often  referred  to  as  the  simplified  motion  equations. 

Specifically  it  is  assumed  that  (1)  the  angle  of  attack  in  pitch 
and  the  roll  angle  are  small  so  that  r  =  i{j  ,  as  in  Ref.  2;  (2)  the 
force  Y,  perpendicular  to  the  torpedo  axis  and  in  the  yaw  plane, 
is  given  by  Y  =  Yp(3  +  Yrr  +  Ygr6r  and  the  moment  N  about  the 
torpedo  center  of  gravity  is  N  =  Npp  +  Nrr  +  Ng^Sj-;  (3)  the  torpedo 
speed  and  mass  are  constant;  (4)  the  angle  of  attack  in  yaw  (P)  is 
small  so  that  sin  p  ar  p  and  cos  p  i.  Eliminating  P  and  substitut¬ 
ing  the  dimensionless  coefficients  (Ref.  2)  leads  to  the  following 
equation  for  the  yaw  angle  vjj; 


...  V  /  Y 
(17)  tii  +  _ 
l 


m  t 


p_ 

J  7.1 


V*l Y'rNp  -  m'LN'p  -  Y'pN'r\  . 


i2 


m-T  J  '7 


V2  NV 


{2  j' 


6v  + 


YpN  6r  -  Y6rN'p| 


m  i 


It  may  be  assumed  without  any  loss  of  generality  that  the  desired 
heading  is  i(j  =  0,  so  that  Eq.  17  is  the  equation  for  the  error  in 
heading.  It  is  then  necessary  to  control  i|j  to  zero.  Assuming 
the  control 
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K 


(18) 


6r  = 


Xc 

TY'«r 


and  referring  to  the  equations 

(19)  Y'p  =  Y’pB  +  L’pT 


Y'r  =  Y‘rB  -  -f-  L*Pt 


N 


Xc 

6r  =  —  Y6r 
r  f 


x. 


NP  “  nPb  +  t  lPt 


M  1  -  TvT 1  /  P  1  T  ' 

N  r  -  N  rg  -  I  I  L 

Eq.  17  becomes 

(20)  ^  +  (a  +  bLpT)i]i  +  (c  +  dLpT  +  eK)4;  +  K(f  +  gL'pT)i|j  = 

where 

v/y'Pb  n'ib\ 


a  = 


i  \  m'-p  j'z 


v 


b  = 


c  = 


i  i 


m  -p 


/t\‘ 
'-*p/  ~i 

J' 


V2  1 


I  _  I 


i  2  m'ipj'.g 


(YrBNpB  “  mLNpB  YpBNr]3) 


d  = 


V  x0  1  x 

“  vl  -Tl  I  ,  P  vl  -Tl 

—  ~ r-T  YrB  -  NpB  -  mL+— YpB  -  — Nr 
i  i  mTJ  z\  *  xp 
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Appendix  B 

DEVELOPMENT  OF  STABILITY  CRITERIA 

The  necessary  and  sufficient  condition  for  stability  of  the 
solution  i}»  =  0  of  Eq.  21,  Appendix  A,  is  that  the  real  part  of  the 
roots  of  m2  +  (a  +  bLp^Jm^  +  (c  +  dLp^,  +  eK)m  +  Kf  =  0  be 
negative.  To  verify  this,  let  the  roots  be  mj,  nr^,  and  m3  and 
consider  the  following  cases: 

1.  If  mp  m2,  and  m3  are  all  real  and  distinct,  then 

,  _  A  mit  .  D  mot  ,  „  mot 
ijj  =  Ae  1  +  Be  c  +  Ce  J 

2.  If  m2  =  m3  and  mj  /  m2,  then 

i|i  =  Aem +  (B  +  Ct)em2t 

3.  If  mj  =  m2  =  m3,  then 

4>  =  (A  +  Bt  +  Ct2)emit 

4.  If  m2  and  m3  are  conjugate  imaginaries  y  ±  Xi ,  then 

4*  =  Aem^  +  (B  cos  \t  +  C  sin  \t)e"'|,t 

In  each  of  these  cases  the  exponent  must  be  negative  in  order  that 
i|j  -*■  0  as  t  -*■  00,  and  conversely,  which  is  simply  the  statement 
that  the  real  part  of  the  roots  is  negative.  By  noting  the  rela¬ 
tionship  between  the  roots  and  the  coefficients,  it  is  easy  to 
derive  equivalent  necessary  and  sufficient  conditions  in  terms  of 
the  coefficients,  i.  e.  , 
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b<jj*  -  2d£,io^  +  (da  -  bc)co^ 

(26)  K  = - - - 

becj“  -  2bf£u  +  df 

and 

-fK 

(27)  m  l  =  — - 


The  author  of  Motion  Equations  for  Torpedoes  (Ref.  3) 
suggested  that  a  response  generally  suitable  for  well- controlled 
torpedoes  is  given  when  4  =  cos  45  0  =  0.707  and  w  =  2tt  .  For 
these  values  let  the  solutions  of  Eq.  24  be  L  p-p  =  L>  *,  K  =  K* 
and  mj  =  m*  .  Then  from  Eq.  25,  26,  and  27 

f(8.886a  -  c  -  39.48)  +  e(350.7  -  39.48a) 

(28)  L*  = - 

f(d  -  8.886b)  +  39.48be 


(29) 

and 

(30) 


K* 


b(  1558  -  39.48c)  +  d(39.48a  -  350.7) 
f(d  -  8.886b)  +  39.48be 


-fK* 


39.48 


The  solution  of  the  yaw- error  equation  for  these  values  of  Lp^, 
K,  and  mj  is 

(31)  \\>  =  A  exp(m*1t)  I  (B  cos  4.443t  I  C  sin4.443t)  exp(-4.443t) 

where  A,  B,  and  C  are  constants  determined  by  the  initial  con¬ 
ditions.  The  term  A  exp(m>;jt)  will  rapidly  approach  zero  if 
-m*j  is  sufficiently  large.  For  instance,  it  will  be  reduced  to 
one  tenth  of  its  initial  value  in  0.2  second  if  -m^  =  11.5.  Hence, 
if  -m^  is  approximately  of  this  magnitude  or  larger,  its  effect 
is  negligible  after  a  very  short  period  of  time,  and  the  torpedo 
will  be  comparable  to  those  considered  in  this  report. 
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Thus  a  torpedo,  designed  to  the  criterion  Lp,p  =  L*  and 
having  the  control  6r  =  [K*/(xg/i ) Yg  ]i{i  ,  will  control  a  yaw 
error  to  zero  according  to  Eq.  31,  provided  the  rudder  is  free  to 
deflect  to  the  angle  prescribed  by  the  control  equation.  However, 
the  rudder  is  often  limited  to  a  certain  maximum  deflection,  say 
±6rQ  (the  "±"  indicates  equal  rudder-deflection  limits  in  both 
directions),  in  which  case  the  torpedo  will  respond  according  to 
Eq.  3  1  if  and  only  if 


(32) 


x 


6  v' 

-iso. 

,  °r  ro 


x 


6  y1  e 

—  Y6jr 

l 


r  1  o 


<  ij/  < 


K* 


K* 


This  may  be  referred  to,  quite  descriptively,  as  the  "band  of 
proportional  control."  When  the  angular  error  falls  outside  this 
band,  the  equation  of  motion  becomes,  for  Lp,^  =  L*  (and  as¬ 
suming  the  error  is  positive,  so  that  the  rudder  is  positioned  at 


Xg 

(33)  4i  +  (a  +  bL*)i|/  +  (c  +  dL*)ib  =  -f  _  y'c  6,. 

ro 


The  necessary  and  sufficient  condition  for  stability  of  the  steady- 
state  solution  of  Eq.  33 


(34) 


c  +  dL* 


is  that  c  +  dL*  is  positive.  For  each  of  the  bodies  considered 
in  this  report  L*  was  sufficiently  large  to  satisfy  this  condition 


46 


NAVORD  REPORT  6428 


and  therefore  this  restriction  is  not  expected  to  affect  the  use 
of  the  criteria. 

Often  the  steady- state  turn  rate  per  degree  of  rudder  deflec¬ 
tion  (^ss/^r)  is  a.  design  specification.  When  6r  ~  Sr0>  from 
Eq.  34,  Yg  is  determined  as 

,  'J'ss 

(35)  Y6  = - 

6r 

o 

and  this  in  turn  will  determine  the  band  of  proportional  control 
from  Expression  32. 


c  +  dL* 


-f 


XS 
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Appendix  C 

SAMPLE  COMPUTATION  OF  CRITERIA 

According  to  Ref.  1  and  2  and  the  report  of  footnote  2,  and  the 
assumptions  of  the  section  "Parameters  of  Sample  Torpedoes," 
Body  1.1  has  the  following  dimensionless  coefficients  and  physi¬ 
cal  characteristics: 

Y'pB  =  0.54,  Y'rB  =  0.30,  NpB=1.10,  N'rB  =  -0.15 

Xp/f  =  -0.446,  x6/f  =  -0.468,  V/t  =  8.906/sec,  j'z  =  0.207 
ztil  =  2.006,  m'-p  =  3.454,  d  =  1.012  ft 

i=  7.590  ft,  in  =  0.304  ft,  lc=  4.554  ft,  it  =  2.732  ft 

Therefore 

ln  0.304 

T|  =  -  =  -  =0.1 

fn  +  lt  0.304  +  2.732 

ip  *n  +  it  0.304  +  2.732 
d  d  1.012 

and 

tc  4.554 

—  = - =  0.6 

i  7.590 

The  quantities  a,  b,  c,  d,  e,  and  f  necessary  to  determine  L*, 
K*.  and  Y'^/^g s/6  rQ)  are  given  by  Appendix  A  as: 


PRECEDING  PAGE  BLANK -NOT  FILMED 
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/  0.54  0.15  \ 

a  =  8.906  - + - I  -  7. 845 

\  3.454  0.207/ 


b  = 


8.906 


- f- 

3.454 


(0.446)2 

0.207 


11.136 


(8.906)' 


c  = 


d  = 


3.454(0.207) 

(8.906)2  (-0.446) 
3.454(0.207) 

(8.906)2 


-[0.30(1.1)  -  2.006(1.1)  +  0.54(0.15)]  =  -199.24 


0.15 


0.30  -  1.10  -  2.006  -  (0.446)(0.54)  - 


0.446] 


=  167.38 


e  =  -- 


=  -383.17 


0.207 


-(8.906)3  /  1.1 

f  = -  0.54  + - 

3.454(0.207)  \  0.468/ 


-2855.3 


Substituting  in  Eq.  7 

-2855. 3[8. 886(7. 845)+  199.24-  39.48]-  383. 17[350. 7  -  39.48(7.845)] 

L*  = - =  1.84 

-2855.3(167.38  -  8.886(11.136)]+  39 . 48(  1 1 . 13 6)( - 3 83 . 1 7) 


Equation  8  gives 

11. 136[  1558  +  39.48(199.24)]  +  167. 38[39. 48(7. 845)  -  350.7] 

K*  = - =  -0.270 

-2855. 3[  167.38  -  8.886(11.136)]  +  39.48(11. 136)(-383. 17) 


From  Eq.  13 

Yg  -199.24  +  167.38(1.84) 

- L_  =  - =  -0.0819 

^ss/5r0  -(2855. 3)(0. 468) 
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If  the  approximate  formulae  are  used,  and,  since  q  =  0.1, 
ip/d  =  3,  and  !c/i  =  0.6,  Eq.  14,  15,  and  16  give 

L*  =  [0.909  -  0.37(0.1)]  (0.6)  +  0.1  +  1.205  =  1.83 

-K*  0.019(3)  -  0.016  +  [0.1  (3)  +  0.082](0.6)  =  0.270 

and 

- - —  *  0.007(3)  +  [0.0185(3)  +  0.0448](0.6)  =  0.0812 

'I'ssAro 
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Appendix  D 

DIMENSIONLESS  TAIL-LIFT  COEFFICIENT  AS  A 
FUNCTION  OF  FIN  SIZE 


The  quantity  Lp,p  is  defined  as  tail  lift/ (l/2)pAV^  and  is 
referred  to  as  the  dimensionless  tail-lift  coefficient.  Tail  lift 
may  be  generated  in  many  ways  (by  shroud  rings,  propellers, 
pumpjets,  etc.  ),  but  the  most  frequent  means  is  cruciform  tail 
fins. 

The  following  drawing  defines  the  necessary  dimensions  for 
relating  fin  size  to  tail  lift: 


a  =  average  body  radius  plus  the  average  momentum  thick¬ 
ness  of  the  boundary  layer  (at  the  fins) 
b  =  the  half- span 
c  =  average  chord 
d  =  maximum  diameter 
From  pages  17  and  46  of  Ref.  2 


2bc  CL 


P 


A 


8bc  CLp 


ltd 


2 


PRECEDING  PAGE  BLANK -NOT  FT )  Aren 
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where 


2trb/  c 


+  1+1 


Therefore,  upon  substitution 


This  relationship  is  approximate  and  means  of  improving  it  are 
being  studied,  but  it  is  accepted  as  reasonably  accurate  for 
preliminary-design  purposes. 
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